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Abstract. In this paper we test the astrometric precision of VLT/FORS2 observations using a serie of CCD frames taken in 
Galactic bulge area. A special reduction method based on symmetrization of reference fields was used to reduce the atmospheric 
image motion. Positional precision of unsaturated R = 16 mag star images at 17 sec exposure and 0.55" seeing was found to be 
equal to 300 pas. The total error of observations was decomposed into components. It was shown that astrometric error depends 
largely on the photon centroiding error of the target (250 pas for 16 mag stars) while the image motion is much less (110 pas). 
At galactic latitudes to about 20°, precision for a serie of frames with a 10 min total exposure is estimated to be 30-50 pas for 
14-16 mag stars providing the images are not overexposed and the filter /? spcc j a i is used. Error estimates for fields with smaller 
sky star density are given. We conclude that astrometric observations with large telescopes, under optimal reduction, are never 
atmospheric limited. The bias caused by differential chromatic refraction and residual chromatism of LADC is considered and 
expressions valid for correcting color effects in the measured positions are given. 
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1. Introduction 

The studies of exoplanets and microlensing effects are usu¬ 
ally restricted to measuring two physical quantities: radial ve¬ 
locities and/or brightness of stars. Angular measurements are 
rarely considered due to their poor precision. The extraction 
of many parameters related to the studied object, from a sin¬ 
gle physical quantity, however, is not always possible and re¬ 
liable. So, radial velocity data allow to find only the low limit 
of planet mass. The additional information on the star angular 
displacements essentially improves an accuracy of exoplanet 
mass and orbital parameters determination (Pravdo & Shaklan 
1 996 Benedict et al. 120021 . This data are also useful for mi¬ 
crolensing studies of massive compact objects in the galaxy 
because allow to measure lens mass, its distance and veloc¬ 
ity (e.g. Boden et al. [1998 1. Safizadeh et al. ( fT999l l argue that 
a combined use of photometric microlensing and astrometric 
measurements allows detection of low mass planets, determi¬ 
nation of their masses and semi-major axes. 

To be of practical importance, an accuracy of astrometric 
measurements for the discussed purposes should be at least of 
the order of 10 to 100 pas. A long history of ground-based 
observations testifies however that a real accuracy is much 
worse, about 1 mas per a night for a 10-20' reference frame 
(Gatewood 1987 ) and 1 mas/h for a 1' double star separation 
measurements (Han T 9E91 l. The best 150 pas/h precision was 
achieved by Pravdo & Shaklan (31 9961 ) in a special serie of ob¬ 
servations at a 5-m Palomar telescope. A major factor limit¬ 
ing the accuracy is the image motion caused by atmospheric 
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turbulence and displayed as a random relative change of star 
image positions in unpredicted directions and at uncertain an¬ 
gle. Therefore development of ultra high-precision astrometric 
methods of observations from the ground is normally related to 
infrared interferometers which have 10 pas expected precision 
(e.g. Frink et al. 2004 ). A similar accuracy, at acceptable expo¬ 
sure times, was found to be unattainable for filled one-aperture 
telescopes due to atmospheric image motion (Lindegren ll980l l. 

This conclusion was revised in the previous work 
(Lazorenko & Lazorenko l2004l further Paper I) where we con¬ 
sidered the image motion as a turbulent light phase-related 
quantity dependent on initial wave-front fluctuations at the tele¬ 
scope entrance pupil, on the process of differential measure¬ 
ments and on the way how these phase fluctuations affect the 
measured positions. In the spectral domain, transformation of 
phase fluctuations into image motion is described by four filters 
of non-atmospheric origin which correspond to: 1) conversion 
from the phase to the wave-front gradient, 2) averaging over the 
entrance pupil, 3) averaging due to a finite exposure and 4) for¬ 
mation of function differences in directions to reference stars. 
Favorable for the inhibition of atmospheric image motion spec¬ 
trum is that circumstance that a filter Y(q) correspondent to the 
averaging over the entrance pupil is a low-pass filter transpar¬ 
ent to about q ~ 2/D spatial frequency (I) is the telescope di¬ 
ameter) while reference field stars form a high-pass filter Q{q). 
The filter bands are only partially overlapping by their rising 
(for Q(q)) and descending (for Y(q)) branches. The shape of 
these instrumental functions can be adjusted so as to minimize 
the combined system response Y(q)Q(q). Improvement of the 
Y(q) filter shape is achieved by apodization, or applying to the 
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entrance pupil a special covering with a variable light trans¬ 
mission. This measure results in a fast attenuation of the filter 
transmission at q > 2/Z). The function Q(q ) shape is adjusted at 
the reduction phase by setting special weights to the each ref¬ 
erence star which virtually reconfigures the reference field into 
a highly symmetric star group and supresses Q(q) responce at 
low frequencies. 

Attenuation factor of the system with modified filters Y(q) 
and Q(q) depends on the telescope diameter, angular reference 
field radius R , the turbulent layer height h and, being a power 
function of the ratio Rh/D is especially large for narrow-field 
mode of observations when Rh <s D. For that reason, the 
method is recommended for future large 30-100 m telescopes 
though for existing 8-10 m instruments the expected precision 
is also high. Based on theoretical asumptions, we had shown 
that for 10-m telescopes precision of differential referencing to 
the background stars is, depending on the brightness and num¬ 
ber of reference stars, from 10 to 60 /ras at 10 min exposure. 

In the current study we test the effeciency of this astromet¬ 
ric method application to actual observations at VLT and in 
the presence of numerous noise sources, such as atmospheric 
image motion, photon noise in the images of stars, pixeliza- 
tion, complex shape of PSF, background noise, optical aberra¬ 
tions, active optics performance etc. In Sect.6 we discuss an¬ 
other atmospheric effect known as differential chromatic re¬ 
fraction (DCR) and which causes a relative displacement of 
star images of different colors (Pravdo & Shaklan ll996t Monet 
et al. fl992l Louarn et al. l2000l . At VLT, the DCR effect is 
strongly reduced by use of a special Longitudial Atmospheric 
Dispersion Compensator (LADC) which mimics atmospheric 
refraction but introduces it in the opposite direction (Avila et 
al. fl997ll . 

2. Observational data and computation of 
centroids 

For the test study we used a 4-hour serie of FORS2 frames 
obtained by Moutou et al. (' 2004 i when observing OGLE-TR- 
132b object with the filter R S peciab at 17 sec average exposure 
and a good FWHM=0.55" seeing. The HR mode with 2x2 
pixel binning and 0.125" pixel size was used. Only master chip 
of CCD with a field of view of 3.9 x 2.1' containing about 400 
stars brighter than 20 mag was considered, and only well ex¬ 
posed stars fainter than 16 mag(R) were measured. Four frames 
obtained with V filter provided color information used for ex¬ 
amination of chromatic effects (Sect.6). 

A full profile fitting was used for computation of x, y co¬ 
ordinates of stellar centroids. The shape of PSF was fitted by a 
12 parameter model of a sum of three elliptic Gaussians with a 
common x, y centre: 

PS F(x,y) - IG(x,y) + (x - x) 2 FG'(x,y)+ 

(y-y) 2 I"G"(x,y) " (1) 

where x, y are pixel coordinates. The main Gaussian G with 
5 free parameters (x, y, gaussian width parameters ctq x , cr o v , a 
term ao specifying orientation of semi-axes) and a flux I (sixth 
parameter) contains about 98-99% of the total star flux; two 


auxiliary Gaussians G' and G" with fluxes /' and I" are ori¬ 
ented along x, y axes. As model parameters, they both include 
x, y, two width terms cr lx , cr lv for G' and similar cro*, <X 2 V terms 
for G". Note that since image motion affects photocenter posi¬ 
tions, it is useful to treat x, y in Eq. 0 not as the ’’fitting model 
center” but as the weighted photocenter defined by equations 
f f(x — x)PSF(x,y)dxdy = f f(y - y)PSF(x,y) dxdy = 0. 
This interpretation emphasizes a necessity to use the PSF mod¬ 
els with first derivates on model parameters (except x, y) being 
symmetric functions of coordinates x, y. In this case small er¬ 
rors in the parameter determination induce symmetric bias of 
the PSF shape which does not affect computation of x, y. For 
that reason, the model 0 does not contain odd coordinate pow¬ 
ers. 

Transition from the continuous expression for the number 
of electrons 0 to the discrete pixel counts requires numeri¬ 
cal integration of PSF(x,y) within the pixel limits. To avoid 
this procedure, we used analytical approximation of the inte¬ 
gral precise to 1CT 4 . 

Solution of Eq. 0 was obtained in square 10x10 px win¬ 
dows which contained main light signal and are small enough 
to minimize overlapping of nearby images. Approximation of 
the actual PSF shape by the model 0 typically is accurate to 
5-10%, the photon noise statistics essentially differs from the 
Poissonian. Deviations of the actual PSF shape from the model 
have a specific wavelike pattern which is highly correlated for 
all stars images of a certain frame (especially at short spatial 
scales) but completely changes at the next frame. At a high 
level of a light signal (bright images, central pixels), these de¬ 
viations strongly exceed random photon noise. 

Using initial fitting for each star, we averaged the ’’mea¬ 
sured PSF - model” residuals over all star images in a frame to 
find systematic part of these residuals as a function var PS F{x— 
x,y - y). To obtain this function estimates at subpixel level, it 
was approximated by a set of 2D cubic polinomials defined in 
3x3 pixel areas with centres displaced at 1 pixel. The func¬ 
tion var/ J .S' Fix — x, y - y) was then subtracted from individ¬ 
ual star images and a final profile fitting applied. This slightly 
improved the precision of the centroiding. One may note that 
the measured PSFs usually show variations at spatial scales of 
some hundreds of pixels. Therefore it seems reasonable to ob¬ 
tain more precise var PS Fix - x, y - y) function shapes by per¬ 
forming averaging of the ’’measured PSF - model” residuals 
only in the limited areas of a frame. Leading to really better 
fitting of PSFs, this consideration is however incorrect with re¬ 
gards to determination of photocenters because application of 
variable image shape corrections dependent on position of a 
star in the frame induces a space-variable bias in the weighted 
center positions. 

An attempt to assign weights to pixels depending on these 
residual values failed also, though the quality of PSF approxi¬ 
mation according to^ 2 criterion improved. This criterion seems 
to be improper for precise centroiding works since inspite of 
low;y 2 ~ 10 fitting quality for bright stars, the actual centroid¬ 
ing precision is near the photon limit (Section 4). 

To estimate the centroiding precision for images with PSFs 
of complex shape 0 , we created random images adding 
Poisson noise to the model function. At 340 e /px sky level 
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typical for FORS2, the results are represented by the expres¬ 
sion 


e = 


FWHM 

2.26 V7 


(1 +475r 0 ' 7 ) 


( 2 ) 


Here FWHM refers to x or y axis and I is given in electrons. 
Note that Eq. 0 is valid for random photon fluctuations only 
and does not take into account a real photon statistics in central 
pixels. It was found (Sect.4.1) that for bright images a limit 0 
is not achievable. 


3. Reduction model 


Astrometric reductions are based on theoretical considerations 
of Paper I which main results are reproduced here in a concise 
form. 

The variance tr 2 t of differential atmospheric image dis- 
placemets caused by a thin turbulent layer can be found by 
integration of power spectrum of image motion which, as it 
was mentioned in Introduction, contains filter functions Y(q) 
and Q(q). Therefore, using Eq.(12) of Paper I, we obtain 


cr: 


m r 

~ VT J 0 


Y(q) 


2^q 2m F 2m (xi,yi) 


n— 1 


q dq 


(3) 


where T is the exposure, V is the velocity of the turbulent layer 
and J(h) is a function proportional to the intensity of a tur¬ 
bulence in the layer and its thickness. The integration is per¬ 
formed in the area of a 2D spatial frequency q related to the 
turbulent layer plane. Cartesian star coordinates x,-,y, not dis¬ 
torted by the atmosphere, refer to the same plane. The subinte¬ 
gral component in square brackets represents the Q(q) function 
and describes filtration caused by reference star field. In origi¬ 
nal form (Eq.(25) of Paper I), it is written as a linear combina¬ 
tion of Bessel zero-order functions dependent on the star coor¬ 
dinate differences. Note a modal structure of the Q(q) compo¬ 
nent expansion into q 2 '" powers with coefficients F 2m which are 
functions of 2 m order cross-moments of star coordinates x,, y, 
differences measured with reference to the target. The ampli¬ 
tude of q 2 m F 2m components rapidly decreases with m. The Y(q) 
function describes filtration of the original image motion spec¬ 
trum by the telescope entrance pupil which for better inhibition 
of high-frequency spectral components can be apodized; for a 
filled pupil an expression for Y(q) is given by Eq.(13) of Paper 
I. 

It is easy to formulate conditions at which the principal 
component F 2 and possibly some next F 2m functions turn to 
zero at a reduction stage. For that purpose, it is necessary to 
define the target position with reference to N field stars by the 
quantity (Eq.(19) of Paper I) 


W = AT 1 


^ cii(x o - Xi) = x 0 - N 1 y atXi 

i= 1 f=l 


(4) 


for the x-axis and by a similar quantity for the y-axis. Here x,, 
y, are the measured (image motion included) Cartesian coordi¬ 
nates of reference stars i = 1,2 ... A; index i = 0 refers to the 
target. Weights — oo < a, < oo are determined from conditions 

E di(Xj - x 0 )“(y,- - yof = 0, a +/3 = 1... § - 1, 

E cij = N, E a 2 sj = min 


where a and f3 are non-negative integers, s,- is the centroiding 
accuracy for i-th star. The first equation in 0 describes for¬ 
mally a system of points with zero weighted cross-moments of 
the k/2 - 1 order and the symmetry center at xo, yo- The even 
parameter k = 2,4... is called therefore the reference field 
symmetry order and is related to the power of the first q k F k 
non-zero component of Eq. 0 . The last condition in 0 is use¬ 
ful for minimization of reference field centroid errors which 
affects the accuracy of W positions. The system 0 is solved 
when at least N m i n = k(k + 2)/8 reference stars are available. 

In the simplest case of k — 2 the first condition of the 
system 0 is not used and weights are found from equations 
E cq — N, Y, fl2 s 2 = min. This refers, for example, to the dou¬ 
ble star (N - 1) separation W x — (xo — xi) measurements or 
to the average of individual separations W x = A' 1 E(xo ~ 4 1 ) 
if N > 1 reference stars of equal brightness are available. The 
field is not symmetrized, and cr~ t value is maximum since the 
largest q 2 F 2 term in Eq. 0 is not eliminated. 

The first k = 4 symmetrized solution is obtained applying 
conditions E a,(x, - xo) = EC V; - Yo) = 0. In particular, these 
conditions refer to the continuous distribution of reference stars 
in a circular field which, as shown by Lindegren (119801 yields a 
strong inhibition of the image motion. The same k = 4 symme¬ 
try order with elimination of the q 2 F 2 term can be implemented 
with use of only 3 arbitrary placed stars. The principal term in 
Eq.0 is proportional to q 4 F 4. 

The highest k — 12 order implementation used in this study 
requires at least (V m j n =21 reference stars to apply conditions 
E a,(xi-x 0 ) = E «;(}'(-}’o) = E at(xj-xo ) 2 = 2 a,(x,-x 0 )(y,- 
yo) = ... = E a i(Xi ~ xo)(y, - y 0 ) 4 = E a,(y, - yo) 5 = 0. This 
results in elimination of F 2 , F 4 , F§, and F 10 terms thus Q(q) 
is proportional to q n F\ 2 . 

In proper motion works, the displacement of the target 
star is found as a difference of W(Ep2 ) - W(Ep 1) at epochs 
Ep2 and Epl. Weights a, are determined from conditions 0 
in which the measured positions at Epl (standard frame) are 
used as unknown true coordinates x, = x,(Ep 1), y ; = yj(Epl). 
This substitution, as seen from Eqs. 0 . leads to the identity 
W(Ep 1) = 0 that defines a conditional zero-point of W values. 
Positions W(Ep 1) define also a zero-point of proper motions 
which are computed as 

Px = W(Ep2) = E a i(xo ~ xd = x 0 - x 0 + 4 E a I {x l - x,) (6) 


Symmetrization of reference fields not only suppress at¬ 
mospheric image motion, but also reduces geometrical field 
distortions, in this respect being equivalent to usual reduction 
models. Let’s consider astrometric reduction model based on a 
polynomial expansion 

Xi - Xi = A 0 + A[’°(x; - x 0 ) + Aj’^y; - y 0 ) + ... = 

k/ 2-1 m a 

Ao + E E Kf(xt - x 0 )“(y,- - yof ( 7 ) 

m— 1 a,/3= 0 , 
a+p—m 


of measured x,, y; coordinates at Ep2 over cross-moments of 
coordinate differences x, -Xo, y, -y 0 at Epl of a standard frame, 
and where the index i = 0 again refers to the target. Equation 
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for y has a similar structure. A residual of conditional equations 
at i — 0 is treated as a proper motion 

Hx = x 0 (Ep2) - xo(Epl)-A 0 (8) 

Note that all coefficients A "f with m > 0 had vanished since 
Eqs. 0 are written for a chosen i — 0 star. This result is com¬ 
pletely equivalent to Eq. 0 since, multiplying Eq. 0 by a, and 
than summing up over i, in view of conditions 0 we find that 
Ao = N X a i(xi - Xi). Thus, symmetrization allows both to 
reduce the image motion and exclude all geometric distortions 
(their change between Epl and Ep2) up to the polynomial de¬ 
gree of a + [i = k/2 - 1 inclusive. And on the contrary, if the 
reduction model 0 includes all modes m < k /2 - 1 with no 
omissions, the least squares solution of systems 0 and 00 
are equivalent. It means that though the use of the model 0 
is aimed at correction of geometric distortions, it reduces the 
image motion equally well. For example, a 6 parameter linear 
model normally used in proper motion works (e.g. Pravdo & 
Shaklan fl996l in fact performs the k = 4 order symmetriza¬ 
tion. 

4. Differential positions 

A one-hour serie of j = 1,2 .. .72 frames obtained near the 
meridian was used for a detailed analysis. A frame j — 1 ob¬ 
tained near meridian was chosen as standard ’’epoch 1” frame. 
In this frame, asuming some star as a target i = 0 and using 
Eqs. 0 , we computed a set of weights a, with which the tar¬ 
get star positions Wj at each other j frame (’’epoch 2”) were 
calculated. Computations were performed at k = 2,4.. 12 and 
reference field radii R from 0.2' to 1.1', whereupon all pro¬ 
cedure was repeated, considering in turn each other star as a 
target. Since the actual image motion for so short time is zero, 
the displacements computed, of course, represent errors of W x , 
W x determination. Because a noticeable part of these errors are 
caused by centroiding errors of reference stars, we performed 
iterative refinement of positions x,-, y,-. For this purpose we car¬ 
ried out calculation of preliminary W x , W y values with k = 6 
and near optimal R - 0.6'. These displacements taken with 
reference to their average (lT f ), (IT,) over the whole serie of 
j = 1,2.. .72 frames were treated as centroiding errors and 
therefore subtracted from x,, y, coordinates of the correspond¬ 
ing z-th star. While reducing random errors of reference star po¬ 
sitions, these iterations introduce an extra atmospheric depen¬ 
dent and spatially-correlated bias. For this reason only a small 
improvement of the final precision was obtained and only at the 
first iteration; no improvement was found for k = 2 at all. 

4.1. Astrometric error components 

Fig 0 represents astrometric errors as residuals Wj - (IT) of 
each frame star position Wj from its average {W) over all 72 
frames. Distribution of residuals computed with R = 0.8', k = 4 
and k — 12 is seen to be symmetric and dependent on the star 
brightness which is especially contrast at k = 12. Errors of W y 
are slightly larger due to the image elongation along the y-axis. 
Note the difference of point scatter at k — 4 and k = 12 which 
were obtained with use of identical sets of reference stars. 



Fig. 1 . Astrometric deviations Wj - (IT) of each frame star po¬ 
sition Wj from this star average position (W) as a function of 
instrumental magnitude; plots correspond to 0.8' reference field 
radius and reduction parameters k = 4 and k — 12 


The measured variance Aj; of each star positions com¬ 
puted with application of some k order was found as A 2 = 
(Wj - (IT)) 2 . It should be commented that At defines also 
a precision of IT value zero-points which are set so that 
W(Epl) = 0. The W(Epl) values are thus shifted from (IT) 
randomly by about Ak r.m.s. 

The estimates of A 2 plotted as a function of R show a large 
scatter of points (Fig0i) however a clear dependence on k. R 
and, of course, on the star brightness (the plot is built for 16.5 
mag star) is seen. To model this dependence, we decomposed 
the Aj; value as a sum 

Ak - 0* + C + *1 s o (9) 

of atmospheric component 0 which at large number of refer¬ 
ence stars is approximated by the expression (Eqs.(36-38) of 
Paper I) 

(An = B k (R/l') bi (10) 

with parameters B * and /?/, dependent on k and considered here 
as free model parameters; centroiding noise of reference field 

°rf = N' 1 a 2 e 2 (11) 

where e, is given by 0 ; and the third component r \eo which 
represents the centroid precision for the target star. The coef¬ 
ficient C < 1 is introduced to reflect the decrease of oyf com¬ 
ponent produced by iterative refinement of star positions, and 
q > 1 is the ratio of the true centroid error to its photon-limited 
estimate 0 ; it was treated as a free model parameter for each 
0.5 mag bin of magnitudes. Solution of Eq.0 produced Bk and 
Ip values, C ~ 0.7, // ~ 1 for stars fainter than 17 mag and, due 
to PSF profile imperfection, q 2 ~ 1 .2 - 1.3 for brighter images. 
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Fig. 2. Dependence of the total astrometric error and its com¬ 
ponents on R at k — 2,4.. 12; a) the total error Aj; of 16.5 mag 
star observation: measured point estimates (dots) and approxi¬ 
mation ([9j (lines); errors of a double-star separation measure¬ 
ments are shown with a shift: a real curve is a factor 2 upwards; 
b) reference field components: oyf (dashed lines with different 
symbols) and <x at (solid lines); Large cicles mark points of line 
intersections at optimal radii R opt ; c) the cumulative reference 
field contribution A a+r . All plots refer to right ascention 


FigEb shows dependence of components C<x r f (units: 
0.001 px = 1 millipixel, hereafter mpx; 1 mpx = 125 //as (mi- 
croarcsec)) and cr at on R at each parameter k value. The ob¬ 
served dependencies qualitatively match theoretic behavior of 
plots shown in Fig. 15 of Paper I and correspondent to D = 10 m 
and T = 10 min. At narrow fields, the dominant component is 


oy-f which decreases with R following an approximate depen¬ 
dence (Paper I) 

FWHM k 

O-rf = -(12) 

2.36 R yjnn 4 

where n is the total light from a lxl ' reference field, R is given 
in minutes of an arc, k/2 is supposed even and star number 
N » A^ m i n . For small N =* 3 - 5Af m j n , the actual values dm of 
cr rf are systematically a factor 1.5-2 above the estimate dm . 
At yet lower N the curves for cr at are not built. At R - 1' the 
value of <x at considerably exceeds <x r f. 

A twinning of curves in FigElrelated to pairs k — (2,4), (6,8) 
and (10,12) is a specific feature of densely populated fields 
since a high sky star density favors reduction due to the field 
averaging effects. At k not multiple 4 the values of <x at and cr,f 
are therefore close to approximate estimates G3 and (1 1 21 com¬ 
puted with k+2. Just for illustration of the field averaging effect, 
Fig|2k reproduces the dependency of astrometric error for mea¬ 
surements of xo - x; separation between two stars of approxi¬ 
mately equal brightness. The dependance observed implies that 
to obtain a good astrometric precision, a separation R between 
two stars should be less than 0.1'; at larger R the accuracy de¬ 
grades very rapidly. This is an example of a marginal single-star 
reference field but note that W x positions © are composed of 
xq - x; unit differences averaged over the field. 


Table 1. Parameters B ,t [//as] and b k calculated for A and Y and 
scaled from Fig. 15 of Paper I 


k 

Bk b k 

computed in X 

Bk b k 

computed in Y 

B k b k 

from Paper I 

2 

1380 

0.9 

1260 

0.8 

1190 

0.64 

4 

780 

1.0 

1010 

1.0 

806 

1.69 

6 

340 

1.3 

420 

1.2 

186 

1.63 

8 

310 

1.2 

370 

1.1 

161 

1.83 

10 

180 

1.8 

230 

1.4 

84 

1.66 

12 

160 

1.9 

200 

1.6 

76 

1.84 


Table 1 compares current parameter B k and hi estimates for 
VLT and their expectations for Chilean observatories (Cerro 
Tololo, Paranal and San Pedro Martir) at typical atmospheric 
conditions. This data was found by scaling plots from Fig. 15 
of Paper I given for a D — 10 m and T - 10 min, to the cur¬ 
rent values using a relation Bj, ~ (airmass) x T l D 3 ( Eq. 
(37) of Paper I). Results given in the last column of Table 1 
generally agree with observed values except some differences 
commented below. By definition 03, coefficient B, t is equal to 
cr at at R — V. Its measured value is approximately equal to the 
scaled expectation at k — 2 and k = 4 but starting from k = 6 
is a factor 2-3 higher. This may indicate an actual intensity 
of atmospheric turbulence at high altitudes which, according to 
the discussion in Paper I, may easily show such fluctuations. A 
smaller than expected measured value of coefficient b k at k — 4, 
6 and 8 implies a higher value of <x at value at small R; unlikely 
it is of atmospheric origin and probably is due to some effects 
caused by the active optics performance. 
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Let us consider a sum A 2 +r = cr 2 t +C 2 o- 2 . that determines the 
properties of a noise from reference field only. This component 
(Fig|2t) has a typical parabola-like shape with a minimum at 
the point R opt where cr 2 t = C 2 cr 2 . Note a gradual decrease of 
A a+1 value with applying large k and a corresponding increase 
of R^. In our case, the minimum is achieved at k - 10,12 
and R ~ 0.9', where A a+r = 1.59 mpx for right ascension and 
A a+r =1.91 mpx for declination. The total astrometric error Ak 
changes with R similarly and also has a minimum at R - P opt . 

Centroid errors ps for the brightest unsaturated stars of 
16 mag are equal to 1.86 mpx (230 /zarcsec) for right ascension 
and 2.22 mpx (280 /./arcsec ) for declination and are of the order 
of best attainable A a+r . The atmospheric error components, Ak 
value for these stars and R opt for each k are listed in Table 2. 


Table 2. Astrometric error components (mpx) and optimal 
field radii P opt at k = 2,4...; last column contains data av¬ 
eraged over k; centroid error rjso is given for brightest stars 



2 

4 

k 

6 

8 

10 

12 

(k) 

rr a , 

2.1 

1.7 

, right ascension 
1.3 1.3 

1.1 

1.1 


C(Tri 

2.1 

1.7 

1.3 

1.3 

1.1 

1.1 

- 

pe 0 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

1.9 

A a + r 

3.0 

2.4 

1.8 

1.9 

1.6 

1.6 

A ( a + r > = 1.1 

Ak 

3.5 

3.0 

2.5 

2.5 

2.4 

2.4 

A < k ) = 2.2 

c a , 

2.8 

2.0 

y, declination 

1.5 1.6 

1.4 

1.4 


Co- rf 

2.8 

2.0 

1.5 

1.6 

1.4 

1.4 

- 

ps 0 

2.2 

2.2 

2.2 

2.2 

2.2 

2.2 

2.2 

A a + r 

4.0 

2.8 

2.1 

2.2 

2.0 

1.9 

A ( a + r > = 1.4 

Ak 

4.6 

3.6 

3.0 

3.1 

2.9 

2.9 

A < k > = 2.6 

^opt 

0 . 22 ' 

0.28' 

0.5' 

0.55' 

0 . 8 ' 

0.9' 

- 


Note that for any arbitrary k value (or equivalent order 
m - k/2 - 1 of a polynomial model 0>, one can find always 
a field size R - P opt at which <x at = Ctr v \ and where atmo¬ 
spheric image motion represents less than a half of the total 
variance Astrometric measurements thus are never atmo¬ 
spherically limited providing that the reduction is performed 
with P opt . This conclusion is valid also for the sky areas with a 
poor sky star density as discussed in Sect.5. 

4.2. Averaging over modes 

Positions Wt calculated with different k and different reference 
field radii P opt are related to distinctive sets of stars and hence 
have almost uncorrelated components oyf. Atmospheric errors 
cr at computed with different k (modes m in Eq. 0) are also 
weakly correlated since the effective turbulent layer heights h 
depend on k. Thus, while at k - 2 all turbulent layers from 5 
to 20 km heights contributes equally to the variance of W k , the 
use of k = 12 effectively cut off turbulence contribution from 
h < 20 km (Fig. 12 of Paper I). It is natural therefore to average 


positions W k with weights P k = l/A 2 +r calculated based on 
Table 1; 

w (k) = Y, PkWk/ Tj Pk (13) 

k k 

The total weight of the average is £ P* = 4 .OP 12 . However, 
taking into account almost equal size of reference fields for 
pairs k = (2,4), (6,8) and (10,12), the effective weight is rather 
P e ff — 2 .OP 12 . The variance of weighed positions O is 

A^>=A 2 a+I . ) + / 7 2 £ 2 (14) 

with a new reference field component 

A< a+r> ~ A -rlatk=12 (15) 

decreased almost twice (Table 2) in comparison to its value A a+r 
at k - 12. 


20 19 18 17 16 




star flux I, electrons 

Fig. 3. Astrometric error as a function of the target flux I. 
Different dots - Ak at A: = 2,4. ..12; open circles - error A/ k ) 
of positions averaged over modes; solid line - approximation 
(1 1 41 of A(k> asuming rj - 1 and with ey given by Eq.0. Upper 
panel - right ascention, lower - declination 

Fig 0shows errors Ak (at R - R opl ) and A(/ f > as a function of 
the target brightness. These estimates were obtained from the 
convergence of the measured positions W k and W/ k ) at different 
frames and after applying a small correction for color effects 
(Sect.6). The model dependence id plotted at 77 = 1 (solid 
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curve) predicts well the measured estimates till R <17 mag. 
At brighter end the actual errors exceed expected at 77 = 1 due 
to non-Poisson noise in star profiles (rf ~ 1.2). Roughly, the 
astrometric error A^) is equal to so at R > 17 mag and to 1.5eo 
for brighter stars. An error component introduced by incorrect 
tying to the reference frame thus is rather small. 

5. Astrometric error at low sky star density 

The results obtained above for high sky star density frames can 
be used to predict astrometric performance of VLT at other 
sky regions. The most important parameter necessary for this 
purpose is the integral light flux n registered from reference 
stars. Its average value for our OGLE field is equivalent to 
no = 12 x 10 6 e~ from a square minute of arc (60 stars brighter 
than R = 20 mag at 1 x 1' area). It essentially affects <T,f value 
which according to F,n. (ll 2 \ goes as <x r f ~ n 1/2 while cr at does 
not depend on n {TOJ. At n ^ no the curves in Fig. 1 relevant to 
the cr jf component are shifted moving the point R opt to a new 
position R' opt = 9R opt . 

light signal n (1 unit = 10 6 electrons) 


0 2 4 6 8 10 



relative light signal n/n 0 

Fig. 4. An error of referencing A< a+r ) as a function of reference 
field total light n: based on the scaling expression d! 6 l at C 2 = 
0.5 (solid curve), the same at C 2 = 1 (dashed line) and found 
from a direct reduction of a sparsed original sample of stars 
(open circles) 

Using Eqs. G3 and o for scaling cr at and oyf compo¬ 
nents at a new value of R' from a requirement of a bal¬ 
ance cr a t = Corf we find 6 = [(coC 2 no) /(C^n)] 1 T 2 L+ 2 ) where 
C ( 2 = 0.5. Note that Eq. ill21 gives systematically underesti¬ 
mated cr rf value for frames with a small number of reference 
stars N ~ (3 - 5)N mln . In view of this we introduce the param¬ 
eter a> > 1 equal to the ratio of the actual tr r ( value given by 
Eq. O to its approximate estimate O- For the sample ana¬ 
lyzed, co =*■ 1 at n > 0 . 2 no and at low n = 0.1 hq it grows from 
1 to 2.5 when k changes from 4 to 12. For R opt and /?<• relevant 
to some k we obtain R' opi = [(ajC 2 no)/(C^n)]' /( 2 hk+ 2 > R opi hence 
the expected error of referencing 

A a+r = [( W C 2 « 0 )/(C‘«)] W( 2 il+ 2 ) A a+r (16) 

Based on this expression and using parameter values listed in 
Tables 1 and 2, we have calculated the expected values of A a+r 


at n = 0.1 ... 1 .Olio in two variants, with C 2 = 0.5 and with 
C 2 = 1. The resulting error that corresponds to the error 
d 1 51 of the average over modes is shown in Fig0] 

To control our scaling approach, we performed direct com¬ 
putations of A(k) with use of the original but artificially sparsed 
star sample and with no iterative improvement of reference star 
positions ( C 2 = 1). Results (open circles) match well the es¬ 
timates (1 1 6b . Note that at n > 0.2no a good referencing at the 
1.5-2 mpx level is maintained, which is comparable to the 
bright star centroid precision. At n ~ 0.1«o due to the decrease 
of a light signal from reference stars, the accuracy of the refer¬ 
encing sharply degrades. Besides, at highest k = 12 order the 
optimum field size 2R' opt = 3.2 - 3.5' is near the CCD size. 

At very low n the star number in the area of the optimal 
radius R = R' opt may be smaller than /V nlin and so be insufficient 
for implementation of the desired k order. This is extremaly 
unfavorable since forces one to use lower k or, even worse, 
performe reduction with not optimal R. Using above expres¬ 
sions for R' t and ,V mm and asuming C - C 0 we find that a low 
limit of reference star density necessary to implement some k 
symmetry order is « min = n 0 [k(k + 2)/(8z moRl pt )] 1 + 1 ,bk a>~ 1,bk 
Asuming oj x 2-3 and taking bj, and R opt values from Tables 1 
and 2 , we find that approximately /; mm - 0.01 no at any k from 
4 to 12. This integral star density (about 0.6 stars brighter than 
R=20 mag per lxl' area) is typical to the Galactic pole (Allen 
fl973ll . 

Note that the unequality N > N m ; n can be easily restored in 
any case by incorporation of some extra faint stars as reference. 
This does not degrade precision, on the contrary, the total flux 
n will increase only and, according to Eq. <x r | improves. 
We conclude therefore that even at low sky star density the key 
problem of astrometric error improvement is related to the tele¬ 
scope light signal limitation but the atmospheric image motion 
(see Table 3, last column). 

Table 3. Centroid precision and astrometric error (fi as) at 
10 min total exposure and different n/no for reduction versions 
with C 2 = 1.0 and 0.5. A ratio of <x at to the total error (last col¬ 
umn) and reference field contribution A a+r (last line) are given 
also 



centr 

red. with C 1 

= 1 

red. with C 1 

= 0.5 

O’at 

R 

error 


n/no 



n/no 


ratio 

mag 

T)e 

0.1 

0.2 

1.0 

0.1 

0.2 

1.0 

% 

20 

300 

308 

304 

301 

305 

302 

301 

0 

19 

167 

181 

174 

170 

177 

172 

169 

1 

18 

98 

120 

110 

103 

114 

106 

101 

3 

17 

65 

95 

81 

72 

87 

77 

70 

6 

16 

42 

81 

64 

52 

71 

58 

49 

13 

15 

26 

75 

56 

41 

63 

48 

36 

24 

14 

16 

72 

52 

35 

60 

44 

30 

35 

(0) 

(0) 

70 

49 

31 

58 

40 

25 

50 


Table 3 gives the centroid precision and astrometric error 
o for 14 - 20 mag (R) stars expected from a reduction of a 
serie of frames with a total T — 10 min exposure and of 30 - 
50 min (depending on observation mode) duration. For bright 
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14-16 mag stars we adopted rp = 1.3, for 17 mag we used 
r\- = 1.2, for fainter stars r\ - 1. The estimates are given for 
reduction versions correspondent to C 2 — 1 and C 2 = 0.5. In 
the last line, the total reference field contribution A' +r which is 
equal to the limiting astrometric precision of very bright stars 
(s = 0) is given. At high sky star density, this limit is 25 yuas 
per a 10 min exposure. The last column contains a ratio of at¬ 
mospheric image motion to the total error; as discussed above, 
this portion cannot exceed 50% by the definition of the optimal 
field size. 

For brightest stars, thus, the expected accuracy of obser¬ 
vations in Galactic bulge is equal to 30-50 yuas at 10 min ex¬ 
posure, degrading to 50-80 yuas at n ~ 0.1 - 0.2«o (2-4 -10 4 
stars to R = 20 mag per square degree which corresponds to 
galactic latitudes of 20-30° (Allen iToTsV ). For stars fainter than 
R >17 mag, the error does not depend on characteristics of ref¬ 
erence field due to a small light signal from the target. 


6. DCR and residual chromatism of LADC 

The effects of DCR was considered elsewhere. In particular, 
Monet et al. ( 19921 studied restriction of the accuracy of proper 
motions derived from the reduction of frames obtained at dif¬ 
ferent hour angles. In this study we use mathematical approach 
similar to that used by Monet et al. (119921 . However, we have 
to take account of joint influence of an atmosphere and of chro¬ 
matic corrector LADC (Avila et al. 119971 used at VLT to im¬ 
prove the quality of images. 

It is known that the relative displacement S of two 
monochromatic images with wavelengths A\ and Ai_ is propor¬ 
tional to a tangent of zenith distance z'-S = 8 tan z where 5 is the 
coefficient of atmospheric dispersion dependent on a difference 
A\ -/L. The corrector reduces this displacement, incorporating 
a negative equivalent displacement -( 8 -d ) tan z where d is the 
residual error of the corrector. The observed relative displace¬ 
ment of the images in a vertical direction therefore is d tan z. 
When a serie of frames is obtained at VLT, the corrector is pre¬ 
set at some fixed zenith distance zl and this setting does not 
change. Therefore 

S = <5(tanz - tanz L ) + t/tanzL (17) 

If Ai is an effective wavelength of a target star and Ao is effec¬ 
tive wavelength of reference field, the measured displacements 
S x , S y are equal to the S projection at the CCD coordinate sys¬ 
tem axis ( x : oposite to right ascension, y : along declination) 

S x = S siny + SS y = So - S cosy + Syo, (18) 

where y is the angle formed by a direction to zenith and the 
y axis, .S'o is the S value in a meridian and ,S' v o, .S' v o are zero- 
points. 

The values of 8 and d were found from a least squares so¬ 
lution of Eq. {3 for each star. As input, we used positions of 
stars in all 280 frames which were calculated with k = 10 and 
R = 0.8'. The results are shown in Fig^as a function of rela¬ 
tive color v — r — (v — r) where v-r is the mean instrumental 
color index of a local star group used as reference for some 




relative color v-r, mag 

Fig. 5. Atmospheric dispersion coefficient 8 (upper panel) and 
the residual LADC bias d (lower) as a function of relative star 
color v — r — (v — r); the size of circles is proportional to the star 
brightness (from 3 • 10 4 to 1.2 ■ 10 6 e “); solid lines - approxi¬ 
mation (HI 


target star. The distribution of point estimates in plots is well 
approximated by linear functions 

8 = -2.6 + (233 ± 8)[v - r - (v - r)] 
d = —(11.0 + 0.9)[v-r- (v-r)] 

where values are expressed in mpx and correspond to atmo¬ 
spheric pressure 744 mb, temperature 11° and 11% humidity. 
Residual scatter of points in Figgis cr( 8 ) - 11.9 mpx for 8 and 
cr(d) = 1.3 mpx for d and a factor 4-5 exceeds formal errors of 
8 and d value detrmination. Only partially it can be explained 
by the errors of color determination (+0.02 mag). Most likely, 
this scatter is intrinsic and can not be adequately modelled us¬ 
ing only v-r colors. 

To check a reliability of the results obtained, we performed 
a control computation of a chromatic refractive index p using 
the expression 

r l "'“ 

p = 346mas 

given by Pravdo & Shaklan ( 1996 ). The integration is per¬ 
formed within the wavelengths limits A min , A max defined by a 


A~ 2 T(A)I(A) dA/ £ T(A)I(A) dA (20) 
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combined transmission T(A) of the filter and CCD and 1(A) is 
the star flux function of wavelength. This expression is nor¬ 
malized to 0° and pressure of 1013 mb. We calculated p value 
for stars with blackbody temperatures corresponding to GO and 
MO spectral types which V-R colors 0.52 and 1.1 are typical to 
the analyzed sample of stars. The difference of p values com¬ 
puted for these spectra and reduced to a unit interval of colors is 
38 mas/mag and corresponds to the color coefficient in expres¬ 
sion {03 for 1 5. After reduction to the pressure and temperature 
at VLT, we obtained 216 mpx/mag which, in limits of errors, 
agrees to the model (I19> value. 

Use of LADC considerably reduses the bias caused by at¬ 
mospheric color effects. At precise presetting of LADC zl = z, 
the bias S is minimized to S = d tan zl with coefficient d 20 
times smaller in comparison to atmospheric coefficient 6 . 

Image displacement induced by DCR effect can be rather 
considerable. For example, for a red star with v - r — 1.43 and 
local reference field average color v - r = 0.82, the amplitude 
of the displacement at hour angles t ~ 3.5 reaches 60 mpx with 
an average 20 mpx/hour rate of image motion (Fig. | 6 jl. For our 
session of observations, LADC was set at zl = 45° while in a 
meridian z m = 37.3°. 



Fig. 6. Chromatic image motion (counter-clockwise and down¬ 
wards) of a red R = 17.5 mag star with v - r = 1.43 over a 
CCD for 4 hours: points - measured positions; solid curve - 
approximation flS 

DCR effect is a cause of proper motion systematic errors 
proportional to the difference of S x ,S y values related to the 
median moments of frames representing Epl and Ep2. In view 
of this, observations should be carried out at a minimum differ¬ 
ence of hour angles At. Let’s estimate restriction on Ar if star 
colors and coefficients of Eqs. m are known. For observations 
in meridian at zl = z m , Eq. o is reduced to S — d tun^m. 
Considering y ~ t, from Eqs. £EB) one can see that the largest 
error cr(S x ) ~ Atcr(d) tan z m dependent on At is expected in S x . 
At the desired accuracy, say <r(S x ) = 0.1 mpx and <r(d) quoted 
above, a typical restriction on Ar is roughly 20 min irrespective 
of a star color. The precision of model G 2 J, thus, is marginal 
to calibrate near meridian observations obtained with the filter 

^special- 


Restrictions are much stronger for observations at large 
hour angles. In this case the main error is caused by the 6 term 
in Eq. and even with use of colours a good calibration of 
DCR displacementes is possible only at a very strong restric¬ 
tion Ar < 1 - 5 min. 

7. Conclusion 

This study is a first practical test of an astrometric method 
based on symmetrization of reference fields. We have shown 
that, in contrast to the common belief, differential astromet¬ 
ric observations with large telescopes (prividing an adequate 
reduction performed) are not atmospheric limited. At optimal 
reduction, the photon centroid noise from reference field is ex¬ 
actly equal to the atmospheric noise. Considering extra cetroid- 
ing errors from the target star, the photon noise component al¬ 
ways dominates over the image motion (Table 3). 

This study allows us to predicate that astrometric precision 
of observations at VLT is 30-50 pas for stars of 14-16 mag 
at 10 min exposure (a session of observations taking about 
40 minutes). This precision is sufficiently good for determi¬ 
nation of precise parallaxes and proper motions of stars, ex¬ 
amination of microlensing effects and exoplanets. Of course, 
the measured parallaxes and proper motions would be relative. 
Moreover, since each target star position is measured with ref¬ 
erence to its peculiar set of background stars, zero-points of 
parallaxes and proper motions are different for different stars. 
This circumstance, however, is not critical for the purposes 
of exoplanet and (in a minor extent) for microlensing studies 
which require only relative astrometric data. If necessary, zero- 
points can be easily reduced to a common system by iterative 
procedure. 

The most essential restriction for VLT is the saturation of 
bright images. Also, sufficiently large light signal necessary 
for a good referencing can be obtained if observations are re¬ 
stricted to galactic latitudes 20-30°. Atmospheric color effects, 
generally rather large, can be compensated by relevant calibra¬ 
tions based on star colors and applying restrictions on hour an¬ 
gles. Long-term systematic errors were not considered in this 
study and may be a source of additional bias. 

Existing filled aperture 8-10 m telescopes are powerful as¬ 
trometric tools due to effective averaging of atmospheric fluctu¬ 
ations of a phase over the aperture and by providing strong light 
signals. For larger telescopes, in view of dependences e ~ D l , 
cr rf ~ I) 1 and <r at ~ I) 5/2 a further increase of the accuracy is 
expected . Interesting possibilities are allowed by apodization 
of an entrance pupil (Paper I) and, in particular, by use of adap¬ 
tive optics which improves both centroiding precision and re¬ 
moves low-frequency components of image motion spectrum. 
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